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ABSTRACT 

We  detent*  tht  tuns  d^errtu  raapooja  of  a  finite  duk  demode  undtr  nmrtjboo*  of  ccauuoi 
applied  Dux  (chrnnnpoienaometry).  Tho  oao  utuboa  ii  recovered  by  rolving  tht  difrcraititt  cquauoni 
in  tht  arailtr  cylindrical  tnordniele  lytun  with  tht  Ilia  of  nntabla  duconunuout  integrals.  In  addition, 
we  pteacat  tht  remit  let  tht  earn  of  linear  sweep  tmperemevy,  which  will  ponarally  be  mom  useful  lor 
controlled  current  taptnmtntt  with  microaiaciitidM.  ■  ■  z*  ■  ■  ,  /  ■  r~  '  t 


INTRODUCTION 

There  has  been  much  recent  work  on  the  analysis  of  mass  transport  to  finite 
electrode  geometries.  One  reason  for  this  activity  ties  in  the  increased  interest  in  new 
applications  of  microelectrodes  [1].  The  applicability  of  the  preferred  spherical 
microeiectrode  geometry  has  so  far  been  somewhat  restricted  eg.,  to  the  elec¬ 
trodeposition  of  ensembles  [2,3]  or  single  mercury  droplets  [4-6],  the  electrolysis  cl 
dispersions  [7,8],  and  the  dropping  mercury  microeiectrode  [9].  Disk,  band,  and  the 
recently  introduced  ring  microelectrodes  [5,10-12],  are  in  general  more  easiK 
constructed,  but  the  necessary  mathematical  analysis  has  so  far  proved  to  be  rather 
intractable.  The  mathematical  difficulties  are  due  to  the  discontinuities  at  the  edze-. 
of  the  electrodes  (e.g.  constant  concentration  or  flux  over  the  surface  of  the 
electrode,  zero  flux  over  the  adjacent  insulating  surface).  The  diffusion  limited  flux 
becomes  infinite  at  these  discontinuities  ( the  combined  effects  of  the  finite  rates  ■ 
the  surface  reactions  and  of  the  distribution  of  potential  and  concentration  acr,  -• 

0022-0728/88/S03.50  1  1988  Elsevier  Sequoia  SA. 


L 


Copy 

permit 


bl«  »o  DT1C  doos  not 

legible  reproduction 


2 


the  surface  i.e..  the  “  tertiary  current  distribution,”  however,  will  limit  the  rates  at 
the  edges  for  real  systems).  Disk  and  ring  microelectrodes  have  the  advantage  that 


quasi-spherical  diffusion  fields  are  established  at  relatively  short  times;  m  contrast, 
diffusion  to  line  or  band  electrodes  does  not  reach  a  steady  state,  the  flux  varying  as 
1  /In  t.  Spherical  diffusion  fields  at  small  electrodes  give  rise  to  high  rates  of  mass 
transport  to  the  surface  so  that  the  kinetics  of  fast  heterogeneous  reactions  and  of 
fast  reactions  in  solution  can  be  studied  under  steady  state  conditions  (e.g.,  see  the 
previous  analysis,  [13]).  While  it  is  clear  that  fast  reactions  can  always  be  studied  by 
decreasing  the  size  of  the  electrode,  there  are  some  applications  that  are  amenable 
to  study  at  “larger”  microelectrodes.  A  variety  of  analytical  and  simulation  proce¬ 
dures  have  been  used  in  attempts  to  develop  adequate  descriptions  of  the  chro- 
noamperometric  and  chronopoten  hometric  responses  at  disk  and  ring  electrodes 
e.g,  see  111.13-261  In  this  series  of  papers,  we  develop  a  general  approach  to  the 
analysis  of  the  non-steady  state,  and  we  apply  the  method  here  to  the  chronopo- 
tentiometric  case  for  constant  and  uniform  flux  over  the  surface.  The  result  is  an 
exact  expression  that  is  valid  at  all  times,  and  for  any  size  disk,  and  is  therefore 
applicable  to  conventional  eiectroanalytical  experiments,  as  well  as  those  at  micro- 
electrodes.  The  approach  to  the  problem  is  based  on  the  properties  of  discontinuous 
integrals  (see  e-g,  refs.  13/25~ahd  which  we  have  extended  from  the  previous 
analysis  (the  prediction  of  tEe  malt 'transfer  coefficients  for  constant  concentration 
and  constant  flux  conditions  in  the  steady  state)  to  include  time  dependent  mass 
transfer.  In  addition,  we  include  the  results  for  a  linearly  swept  current  experiment. 
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THEORETICAL  CONSIDERATIONS 


General  time  dependent  solution 

For  any  simple  electrochemical  experiment  involving  a  single  reactant,  we  must 
solve  the  time  dependent  diffusion  equation  in  circular  cylindrical  coordinates; 


3c 
3 1 


D  3c 
r  dr 


(1) 


where  c  is  the  concentration  of  the  reactant,  and  r  is  the  radial  distance  coordinate 
measured  from  the  center  of  the  disk  electrode  which  is  imbedded  in  the  insulating 
plane  at  z  —  0.  The  general  initial  condition  is,  at 

r>0,t>0,tm0  C“c"  (2) 

where  c*  is  the  bulk  concentration.  Laplace  transformation  of  eqn.  (2)  gives 


3jc 

3rJ 


+ 


I  —  d:c 

r  3r  + 
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where 


q1  —  s/D  Ml 

and  s  is  the  Laplace  transformation  variable.  The  last  term  on  the  LHS  of  eqn.  (3)  is 


) 


a  constant  and  hence  determines  the  particular  integral  of  eqn.  (1).  We  then  seek  the 
solution  for  the  complementary  function  from 


3:?cf  ^  1  3?cf  +  i2£F  _  x-  „  o 
dr1  r  3r  ds1 

(5) 

Separations  of  vanables  of  the  form 
cCF-uexp(-/(X.  <?);] 

simplify  the  differential  equation  (3)  to  the  familiar  Bessel  form 

(6) 

d2D  1  dti 
dr  2  r  dr 

(7) 

where 

<*J -(/(*•  q)] 2 

For  i, .stance,  we  can  choose  the  simple  form 

(8) 

/(X.  q)  -  (a2  +  q2)i/2 

(9) 

so  that 

cCF-6exp[-(a2  +  q2),/Isj  t. 

(10) 

and  eqn.  (7)  is  recovered  from  eqn.  (3)  with  a  independent  of  q.  (The  factors 
involved  when  considering  a  choice  for  this  function  will  be  discussed  elsewhere.) 
Therefore,  the  solution  to  Bessel’s  differential  equation  (3)  becomes 


?“h m~f  *(*. f)*)4>(«’)  <*«  (n) 

5  JQ 

where  J0  is  the  3~ssd  function  of  the  first  kind,  order  0,  and  we  choose  g( A,  <7)  to 
satisfy  the  boundary  conditions. 

Chronopolentiometry 

We  consider  here  the  solution  for  eqn.  (11)  for  the  chronopotentiometry  problem 
for  a  constant  a  constant  uniform  flu*  -g(mol  cm-1  s'1)  over  the  surface  at  all 
t  >  0.  The  boundary  conditions  at  the  surface  of  the  electrode  become  \ ^  ^ 

0  <  r  <  a,  r  -  0,  f  >  0 

r>  a,  r  —  0,  r  >  0 

We  take  the  Laplace  transform  of  e<jns^£l2),  and  substitute  in  eqn.  (11)  after 
differentiation  under  the  integral  and  obtain,  for  t  >  0. 

ar)  dor  0<r<a,z»0  (1?) 


-  -Q 
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The  boundary  conditions  suggest  that  we  can  apply  the  discontinuous  integrals 

|°  r>a  (15) 

j  J0(ar)J](aa )  da  “  ^  l/2a  r»a  (16) 

°  \  1  /a  r  <  a  (17) 

in  the  solution  of  the  problem.  We  determine  the  conditions  for  which  g(X,  q)  fits 

the  boundary  conditions  (13)  and  (14).  From  eqa.  (13),  we  see  that 

/  g(X,  <j)/(X,  q)J0(ar)  da- 0<r<a,  r-0  (18) 


^jf  *(*.  <?)/(*.  ?)Jo(ar)d«-^ 


Therefore  if  we  let 

,,  x  -Man) 
g(  ’^“/(X.*) 

the  complete  solution  to  eqn.  (3)  having  boundary  conditions  (13)  and  (14)  is 

~c~t-  fsC  exp[-/(x-  «>*i  /o(ar)y'(aa)7(r^ 

and  at  i  -  0, 


.  c*  Qa  r°°  r  /  ^  ,  x  d« 

c’t-dJJq  y«(ar)/>(<w)7(x^) 


We  point  out  here  that  the  interpretation  of  the  arguments  of  the  Bessel  functions 
depends  on  the  nature  of  the  assumption  of  the  form  of  /(X,  q)  and.  indeed,  on  the 
nature  ofjhe  experiment  (e.g.  compare  with  the  discussion  of  the  chronoamperomet- 
ric  case((31j)  The  use  of  the  simple  form  (9)  gives 

-T-* 

for  galvanostatic  conditions  with  a  now  independent  of  q:  we  can  therefore  invert 
immediately  to  the  (-domain: 

Mar)J\(aa)  erf(  Z>l/:arl/2)—  (24) 

where  erf(  y )  denotes  the  error  function. 
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We  can  evaluate  the  average  concentration  at  :  —  0  by  integrating  eqn.  (24)  over 
the  surface  of  the  disk: 

(25) 

u  •'0  a 

With  the  substitutions 


I2  -  Dt 
P  “  a/ 

eqn.  (25)  can  be  written  in  terms  of  dimensionless  variables  and  parameters 


(26) 

(27) 


(28) 


The  function  ♦l  is  tabulated  in  Table  1  as  a  function  of  the  dimensionless 
parameter  (  Dt/ a 1 ). 

At  long  times  and  sufficiently  small  values  of  the  flux,  we  do  not  observe  a 
transition  time  and  always  reach  the  steady  state  value  [13] 


«a*  -  c" 


8ga 

JirD 


(29) 


If  Q  is  sufficiently  large  we  will  get  a  sharp  transition  time  as  the  surface 
concentration  of  the  reactant  approaches  zero;  for  this  condition  eqn.  (28)  can  be 
written 


-  1Qa 


(30) 


from  which  the  transition  time  can  be  obtained.  The  problems  to  be  discussed  in 
this  series  of  papers  are  always  cast  in  the  form  of  such  definite  integrals,  and  these 
frequently  converge  slowly.  Accurate  evaluation  is  readily  obtained,  however,  through 
the  use  of  standard  Bulirsch-Stoer  numerical  integration  methods.  In  addition, 
convergence  can  always  be  speeded  by  suitable  rearrangements.  For  instance,  eqn. 
(30)  can  be  rewritten  in  the  equivalent  form 


2  Qa  If  r»f  r»[  ,/0£\ 

Dc™  a\J0  I'M  l  jj  pi  J0  pi  /  ) 


2(erf(/i)  -  1)~7 

p 


(3D 


The  value  of  the  first  integral  is  known  to  be  4a/3n7  -(2B).  The  second  integral 
converges  rapidly  since  the  error  function  approaches  1  with  increasing  /?.  Figure  1 
illustrates  the  square  root  of  the  dimensionless  transition  times  Dr/a2  as  a  function 
of  the  inverse  of  the  dimensionless  flux  2Q  a/Dcx.  The  values  are  close  to  those 
predicted  by  Aoki  and  Osteryoung  (18]  even  though  these  authors  based  their 
analysis  on  a  uniform  surface  concentration  boundary  condition:  the  expressions 
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TABLE  1 


Value*  of  the  function  <$,(  Dt / a:) 


Dt/al 

Di/a * 

♦i(Pr/a:) 

4.0000  x  10* 

2.1217x10'* 

2.7778 

2.0533x10'’ 

2.7778x10* 

2.5460x10'* 

2.0408 

2.3064  x  10 

2.0408  x  10* 

2.9702x10'* 

1.5625 

15372X10'1 

1.S625  x  10* 

3.3944x  10'* 

1.2346 

2.7473x10-' 

1.2346x10* 

3.8186x10“* 

1.0000 

19381  X  10' ' 

1.0000x10* 

42427x10'* 

2.5000X10'1 

4.1045x10'' 

2.5000x10* 

8.4826X10'* 

UlllxlO'' 

4.5961x10*' 

1.1111  x  10* 

1.2720x10'* 

6.2500X10'* 

4.8524X10'1 

6-2500  x  10* 

1.6954x10'* 

4.0000X10'* 

5.0085x10“' 

4.0000X10* 

2.1185X10'* 

17778X10'* 

5.1132x10-' 

2.7778X  10* 

2.5414x10“* 

10406X10'* 

5.1883x10-' 

2.0406  X 10* 

2.9640x10'* 

1.5625x10'* 

5.2446x10'' 

1.562JX10* 

3.3863x10'* 

13346x10'* 

5.2888x10-' 

1.2346x10* 

3  8083X10'* 

1.0000x10'* 

5.3240X10'1 

1.0000  X 10* 

4.2300x10'* 

15000X10'* 

5  4828  X 10' 1 

2.5000x10* 

8.4318X10'* 

1.1111X10'* 

53358  x  10' 1 

UlllxlO* 

1.2605X10'* 

6.2500X10'* 

5.5623  x  10“ 1 

6.2500x10* 

1.6751X10'* 

4.0000X10'* 

53783X10'* 

4  0000x10* 

2.0868X10'* 

17778X10'* 

5.5889  X 10' ' 

2.7776x10* 

14957x10'* 

10406X30'* 

53964X10-' 

2.0408x10* 

19018x10'* 

13625x1%'* 

5.6021X10-' 

1.5625X10* 

3.3051X10'* 

13346X10'* 

5.6065X10'' 

1.2346x10* 

3.7055X10'* 

1.0000X10'* 

5.6101X10-' 

1.0000x10* 

4.1032X10'* 

15000X10'* 

5.6260X10'* 

2.5000X10* 

7.9259X10'* 

1.1111X10'* 

5.6313X10'1 

UlllxlO1 

1.1472X10'1 

63500X10'* 

5.6340x10-' 

6.2500 

1.4749  X 10' 1 

4.0000X10'* 

5.6356X10-' 

they  derived  could  not  be  inverted  exactly  over  the  entire  time  range  whereat  eqn. 
(30)  is  exact.  It  should  be  noted  that  the  constant  flux  condition  is  more  likely  to 
apply  for  most  of  the  duration  of  the  experiment  rather  than  the  constant  surface 
concentration  condition,  the  actual  behavior  lying  between  these  two  limiting 
conditions.  A  more  correct  formulation  of  the  boundary  condition  on  the  disk  for 
an  irreversible  reaction,  say,  would  be 

8c/3z  —  kc  0<r<a,  r-0  (32) 

For  a  cathodic  reaction,  say,  we  would  write 
k  -  k0  exp[ aEF/RT ]  -  constant  0  <  r  <  a 

provided  we  neglect  the  effects  of  the  distribution  of  potential  in  the  solutio 
we  neglect  the  primary  or  secondary  current  distribution.  Equation  (32)  she 
we  cannot  strictly  speaking  assume  cither  the  flux  or  the  '•oncen'ratio.i  to  K* 
constant  over  the  surface  of  the  disk.  However,  the  close  agreement  of  the  transition 


ik'  ■ e 
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Fig.  1.  Plot  of  the  square  root  of  tht  dimautoolaa  transition  tuna  (Dr/a!)1/5  as  a  function  of  the 
inverse  of  the  dtmeeaonkM  flux  (2Q*/Dcm\ 


times  derived  using  there  two  approaches  shows  that  the  interpretation  is  not  very 
sensitive  to  the  native  of  the  assumptions;  subsequent  paperrfSHiWiU  show  that  this 
is  also  the  case  for  other  types  of  experiment  We  note  mat  the  constant  flux 
condition  will  hold  at  low  current  densities  (as  in  relaxation  experiments)  while  the 
constant  concentration  condition  may  be  approached  in  the  limiting  current  region 
(see,  however,  below)  so  that  the  actual  behavior  must  lie  between  these  two  limits. 
Thus,  while  it  would  be  possible  to  develop/the  conditions  (31)  and  (32),  it  is 
unlikely  that  the  accuracy  of  the  experimenuudata  would  allow  an  assessment  of 
the  range  of  validity  of  the  various  assumptions.  Furthermore,  the  application  of 
eqn.  (31)  would  also  require  the  consideration  of  the  distribution  of  the  potential  in 
the  solution,  i.e.  we  have  to  consider  the  tertiary  current  distribution.  We  can 
predict,  to  some  extent,  the  probable  outcome  of  such  analyses:  it  is  unlikely  that 
the  “throwing  power”  of  any  practicable  system  would  ever  be  sufficiently  low  that 
the  current  density  could  deviate  appreciably  from  uniformity  over  the  surface  of  a 
microdisk  (note  that  an  infinite  flux  to  the  edge  of  the  disk  is  clearly  impossible 
even  on  the  limiting  current  plateau!).  We  therefore  consider  that  the  constant  flux 
boundary  condition  will  hold  under  most  experimental  conditions. 

Figure  1  shows  the  expected  linear  dependence  at  high  values  of  2 fiQa/Dc^ 
where  we  observe  essentially  linear  diffusion  to  the  electrode  followed  by  a  rapid 
rise  at  low  values  of  2 Qa/Dcx  as  this  parameter  approaches  the  condition  required 
for  the  observation  of  a  steady  state 
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The  calculated  values  of  c.„  may  be  used  to  derive  the  potential- time  curves  for 
-rpropnate  models  of  the  electrode  reactions.  For  instance,  in  the  case  of  simple 
Butler-Volmer  kinetics,  we  obtain 


FQ 

'o 


ZQa  .  Di  \  I  -ar\F 


\  .  ( 1  -  a )  ijF )  \ 

rexp(— 


■  exp| 


■avF  i 


xp| 


(1  -  a)yF  | 


RT  I  ^  RT 
which  shows  that  the  transients  are  a  function  of  a,  FQ/i0  and  2 Qa/Dcx . 


(35) 


Linear  sweep  amperometry 

It  is  clear  that  since  there  is  rapid  attainment  of  steady  state  diffusion  to 
microdisk  electrodes,  constant  current  experiments  may  not  be  generally  useful.  It  is 
difficult  to  determine  the  appropriate  galvanostatic  condition  that  will  give  a 
transition  time  within  a  convenient  experimental  time  scale.  It  is  therefore  more 
straightforward  to  apply  time  dependent  fluxes  to  the  surface.  It  we  consider,  for 
example,  the  simplest  case  of  a  linear  current  ramp 


<?(')-  V 

we  can  immediately  rewrite  eqn.  (23)  in  the  general  form  (1,13 
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and  by  the  same  methods  obtain  the  transition  rime 
4  yia 
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or 

4-yra 

~Dc* 
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If  now  the  current  is  swept  from  zero,  a  very  sharp  transition  is  observed  in  the 
potential-tune  plot.  For  simple  Butler-Volmer  kinetics,  the  shape  of  the  response 
can  be  determined  from 
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Table  2  gives  values  of  as  a  function  of  Di/a2. 

The  application  of  other  boundary  conditions  representing  other  electrochemical 
experiments  (e.g.  cyclic  amperometry,  ac  impedance  measurments.  and  to  cases 
involving  reacuons  in  solution  coupled  to  the  electrode  processes)  to  this  form  of 


table  : 

Values  of  the  funcuon  <J MDr/a2) 


[3  <3 


Dr /a2 

Dt/a 2 

4>i<  Dt/a') 

LOOOOx  10~; 

’5112x10-* 

5  0000 

3.6511  x  10'1 

N 

3.1623  x  10' 2 

1  8778  x  10- 4 

7.0711 

3.9321  x  10“' 

] 

3.5355  x  to-2 

2.3473x10-' 

10000 

4  0842x  10“' 

3.7796  x  10  ‘ 2 

2  6826x  10-' 

1.1180 

4  1155x10“' 

j  tC-O. 

4  4721  x  10' 2 

3  ^556x10*' 

1  1952 

4,1313x10“' 

i  viLa  Tit-  ^ 

5  OOOOx  10': 

4  6945x  10*' 

1.4142 

4.1631X10*' 

c 

f 

'J 

7  0711  X10*: 

9  3831  x  10-' 

1.5811 

4.1792x  10*' 

l.OOuOxlO*1 

1  8728x10-* 

2.2361 

4  2116X10*' 

fV  Vt  VYcVatV 

1.1180  x  10" ' 

2.3385x  10-' 

3.1623 

4.2277x10*' 

1.1952x10-' 

2.6706x10** 

3.7796 

4.2327x10*' 

1.4142x10-' 

3.7299x10*' 

4.4721 

4.2360  X 10  *' 

1.5811X10-' 

46526x10*' 

50000 

4.2376X10“' 

2.2361  X  10" 1 

9.2073X10*' 

7.0711 

4.2409X10“' 

3.1623x10-' 

1  8023x10*' 

1.0000X10* 

4.2425  X 10* 1 

3.5355  x  10-' 

2.2285x10“' 

1.1180X10* 

4.2428  X 10  “* 

3.7796x10-' 

2.5269x10“' 

1.1952X10* 

4.2430  X 10“ 1 

4  4721x10-' 

3  4486x10“* 

1.4142X10* 

4.2433X10“' 

5  0000X10'1 

4.2138X10“* 

1.5811X10* 

4.2435  x  10“' 

7  0711x10-' 

7  4851x10“' 

22361x10' 

4.2438  X 10  “' 

1  0000 

1.1751X10“' 

3.1623X10* 

4.2440X10“' 

1.1180 

1.3151x10“' 

3.5355X40* 

42440X10*' 

1.1952 

1.3971x10“' 

3.7796X10* 

4.2440X10*’ 

14142 

1.5820  X 10“ 1 

4.4721  X 10* 

422440X10*' 

1.5811 

1.7040x10“* 

50000X10* 

4.2441X10“' 

2.2361 

22654x10“' 

7.0711X10* 

4.2441x10“' 

3.1623 

19695x10“' 

1.0000x10' 

4  2441X10“' 

3.5355 

3.1727X10*' 

3.1623X10' 

4.2441X10*' 

4.4721 

3.5219X10“' 

1.0000X10* 

42441x10“’ 

analysis,  as  well  as  extension  of  the  analysis  to  include  the  ring  geometry  will  be 
discussed  elsewhere 
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GLOSSARY  OF  SYMBOLS  USED 

a  Radius  of  disk,  cm 

c  Concentration,  mol  cm  “ 3 

cx  Bulk  concentration,  mol  cm'5 

c*v  Average  concentration,  mol  cm-5 

cCF  Concentration  complementary  function 

c'  Surface  concentration,  mol  cm'3 
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D  Diffusion  coefficient,  cm;s_l 

£  Electrode  potential,  V 

£  Faraday  constant.  96485  C  mol*' 

z0  Exchange  current  density.  A  cm  * 1 

J0.  /,  Bessel  functions 
k  Heterogeneous  rate  constant,  cm  s*1 

k°  Heterogeneous  standard  rate  constant.  :m  s*1 
/  <  £>/>'* 

q  (s/D)1'1 

Q  Flux,  mol  cm"2  s"1 

R  Gas  constant.  8.314  J  mol*1  K.*1 

r  Radial  coordinate,  cm 

s  Laplace  transform  variable 

:  Time,  s 

T  Temperature,  K. 

v  Concentration  amplitude 

z  Coordinate  normal  to  plane  of  disk,  cm 

a  Transfer  coefficient  (when  in  exponent) 

a  Continuous  dummy  integration  variable 

0  oU 

y  Flux  sweep  rate,  mol  cm*2  s“l  ’  t 

X  Continuous  dummy  integration  variable 

7j  Overpotential,  V 

*■  jjCNt)] 

r  Transition  time,  s 
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